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Abstract 

We present numerical and analytic results for uniaxial and biaxial order at the isotropic-nematic 
interface within Ginzburg-Landau-de Gennes theory. We study the case where an oblique anchor- 
ing condition is imposed asymptotically on the nematic side of the interface, reproducing results 
of previous work when this condition reduces to planar or homoeotropic anchoring. We construct 
physically motivated and computationally flexible variational profiles for uniaxial and biaxial or- 
der, comparing our variational results to numerical results obtained from a minimization of the 
Ginzburg-Landau-de Gennes free energy. While spatial variations of the scalar uniaxial and biaxial 
order parameters are confined to the neighbourhood of the interface, nematic elasticity requires 
that the director orientation interpolate linearly between either planar or homoeotropic anchoring 
at the location of the interface and the imposed boundary condition at infinity. The selection of pla- 
nar or homoeotropic anchoring at the interface is governed by the sign of the Ginzburg-Landau-de 
Gennes elastic coefficient L2. Our variational calculations are in close agreement with our numerics 
and agree qualitatively with results from density functional theory and molecular simulations. 

PACS numbers: 42.70.Df,67.30.hp,61.30.Dk,61.30.Hn 
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I. INTRODUCTION 



Nematic liquid crystals, typically formed in suspensions of rod-like molecules whose aspect 
ratio deviates sufficiently from unity, exhibit orientational order in the absence of transla- 
tional order Such order is quantified through a traceless, symmetric tensor Q^p 

defined at every point in space l|, Q]. In the nematic phase, the order parameter is 

35* / 1 \ T 

Qafs = ( '^'^'^P " 3 ) + ^ ~ mamp) (1) 

where the director n is defined as the normalized eigenvector corresponding to the largest 
eigenvalue of Q, the subdirector 1 is associated with the sub-leading eigenvalue, and their 
mutual normal m is obtained from n x I. The quantities S and T represent the strength of 
uniaxial and biaxial ordering: |S1 7^ 0, T = is the uniaxial nematic whereas S,T ^ with 
T < 3S defines the biaxial casefll]. 

The description of the early stages of phase-ordering upon quenches from the isotropic 
phase, the properties of nematic droplets within the isotropic phase and the structure of the 
isotropic-nematic interface are all problems which require that nematic and isotropic phases 
be treated within the same framework. The inhomogeneous order parameter configurations 
obtained in these cases are weighted by the Ginzburg-Landau-de Gennes (GLdG) free energy, 
obtained via a gradient expansion in Q in which only low-order symmetry allowed terms are 
retained . The simplest of the problems above is that of the structure of the infinite, 
fiat isotropic-nematic interface, studied initially by de Gennesjsi. 

Nematic ordering is strongly infiuenced by confining walls and surfaces, which impose 
a preferred orientation or "anchoring condition" on the nematic state. Such a preferred 
orientation yields an anchoring angle, defined as the angle made by the director in the 
immediate neighbourhood of the surface with the surface normal. Anchoring normal to the 
surface is termed as homoeotropic, whereas anchoring in the plane of the surface is termed 
as planar. The general case is that of oblique anchoring. 

As is the case with surfaces, the interface between a nematic and its isotropic phase can 
also favour a particular anchoring. The problem of interface structure for the nematic is 
particularly interesting since it illustrates how the structure in the interfacial region can differ 
substantially from structure in the bulk. It is known, for example, that a region proximate 
to the interface can exhibit biaxiality within the LGdG theory, even if the stable nematic 
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phase is pure uniaxialQ], provided planar anchoring is enforced. Such biaxiahty is absent if 
the anchoring is homoeotropic [sl . These two hmits, of homoeotropic and planar anchoring, 
lead to interface profiles of S and T which vary only in the vicinity of the interface, as well 
as orientations which are uniform across the interface js]. 

Can oblique anchoring be stabilized, within GLdG theory, at the interface between a bulk 
uniaxial nematic and its isotropic phase? Suppose we introduce boundary conditions that 
impose a specified oblique orientation deep into the nematic phase, where the magnitude of 
the order parameter is saturated. The question, then, is whether such an imposed orientation 
is relaxed to a preferred value in the vicinity of the interface. The difficulties with this 
problem stem from the fact that changes in the local frame orientation on the nematic side 
of the interface come with an elastic cost arising out of nematic elasticity. This is an effect 
sensitive, in principle, to system dimensions, since gradients can be smoothed out by allowing 
the changes to occur over the system size. While this cost can be reduced by suppressing 
the order parameter amplitudes in regions where order parameter phases vary strongly, the 
precise way in which this might happen, if at all, is an open question. 

Popa-Nita, Sluckin and Wheeler (PSW) ^ studied this problem numerically within a 
GLdG approach, using a set of variables rjs and introduced in Ref. j^. These variables 
are combinations of the variables S, T and 6 used in this paper. Although the focus of their 
study was the emergence of biaxiality at the interface with a planar anchoring condition, 
PSW remarked that if the asymptotic orientation of the director in the nematic phase was 
set to any value other than 90° (planar anchoring) or (homoeotropic anchoring) for large 
z, then rjs and fis approached this value with non-zero slope. PSW thus concluded that 
there could be no stable anchoring if the orientation of the director in the nematic phase 
was neither planar nor homoeotropic, but oblique. The precise nature of the resulting state 



obtained upon applying an oblique anchoring condition was not addressed by PSW[ 

Density functional calculations on hard-rod systems using Onsager's theory applied to the 
free isotropic-nematic interface indicate that the minimum surface free energy is obtained 
when the rods lie parallel to the isotropic-nematic interface, the case of planar anchoring ^, 



10|. Molecular simulations of a system of hard ellipsoids, in which an anchoring energy 



fixes the director orientation in the nematic phase at a variety of angles, indicate that the 
isotropic-nematic interface favours planar anchoring. These simulations, and a mean-field 
calculation based on the Onsager functional, find that the angle profile is approximately 
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linear as one moves away from the boundary condition imposed by the wall at one end of 



the simulation box 
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12| . These results, in particular concerning the stability of planar 
anchoriner, are 

several other papers indicate specific regimes in which homoeotropic or oblique anchoring 
may be stable. Moore and McMullen[l8|] numerically evaluate the inhomogeneous grand 
potential within a specific approximation scheme finding that planar anchoring is preferred 
at the interface for long sphero cylinders, but oblique or homoeotropic anchoring may be an 
energetically favourable alternative for smaller aspect ratios. Holyst and Poniewierski study 
such hard spherocylinders in the Onsager limit, noting that oblique anchoring is favoured 
over a considerable range of aspect ratios {l9|. Finally, experiments provide evidence for both 



oblique |20|] and planar anchoring [2 1|, with electrostatic effects possibly favouring oblique 
anchoring. 

This paper studies the isotropic-nematic interface within GLdG theory in the case where 
an oblique anchoring condition is imposed on the nematic state far from the location of 
the interface. For a flat interface, the components of Q can depend only on the coordinate 
perpendicular to the interface. We assume that this coordinate is aligned along the z axis, 
as shown in Fig. [H which defines the geometry we work with in this paper. We work at 
phase coexistence, imposing boundary conditions fixing the isotropic phase at 2; = —00 and 
the nematic phase at 2 = 00. The components of Q as z ^ 00 are chosen so that 5* is 
fixed to its value at coexistence ^c, while the axis of the nematic is aligned along a specified 
(oblique) direction. The coexisting states must be separated by an interface in which order 
parameters rise from zero on the isotropic side of the interface to saturated, non-zero values 
on the nematic side. Since the two free energy minimum states are degenerate in the bulk, the 
position of the interface is arbitrary and can be fixed, for concreteness, at z = in the infinite 
system. However, there are subtleties. Provided all components of Q vary substantially only 
in the neighbourhood of the interface, the interface can be located through several, largely 
equivalent criteria. However, if variations of Q are not confined to a region proximate to the 
interface but depend on the system size L irrespective of how large L is, the very isolation 
of an interface from the bulk is ill-defined. As indicated earlier, it is this situation which 
obtains in the case of oblique anchoring and the L — 00 limit must be taken with care. 

The central results of this paper are the following: A numerical minimization of the GLdG 
free energy which imposes a specific oblique anchoring condition on the system deep into the 
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FIG. 1: The interface geometry and the coordinate system used in this paper. Note that the 
nematic director makes an angle 9 with respect to the z— axis direction. This angle is fixed at 
infinity io 9 = 9e ■ It can be chosen to vary between 9^ = Q (homoeotropic anchoring at infinity) 
and 9e = 90° (planar anchoring at infinity). The isotropic phase is favoured, through boundary 
conditions, as z ^ —oo, whereas the nematic phase is favoured for z ^ oo. The plane of the 
interface is the x — y plane, shown by ABCD in the figure, whereas the director is confined to the 
EFGH plane as shown. The origin is denoted by O. 

nematic while fixing the interface location at the origin shows that the elements of Q vary 
with space even far away from the interface, albeit slowly. Only in the limit of homoeotropic 
or planar anchoring is the variation of Q confined to a finite region. This variation in the 
case of oblique anchoring can, however, be split into hydrodynamic and non-hydrodynamic 
components. Generically, the variation of the non-hydrodynamic components, such as the 
magnitudes of S and T, are confined to a finite region, independent of the system size L, if 
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L is large enough. However, the orientation of the nematic director varies in space: if the 
asymptotic value of the nematic order parameter at L represents uniaxial ordering along 
an oblique axis, the director orientation interpolates linearly between either a 90° value 
preferred at the location of the interface (planar anchoring) or a 0° value (homoeotropic 
anchoring), and the value imposed by the boundary condition at L. Whether planar or a 
homeoetropic anchoring is preferred at the interface depends on the sign of the second of the 
elastic coefficients in the GLdG expansion, the L2 term, as initially shown by de Gennesjsl. 

Our results are consistent with the qualitative observations of PWS, but provide a detailed 
quantitative analysis in the case of oblique anchoring. We scale angle profiles computed for 
different values of the system size L onto a universal curve, indicating a linear profile. In the 
limit that L — > 00, the slope with which the phase varies vanishes as 1/L, so that the total 
energy cost for elastic distortions of the nematic field ~ J (VOydz ~ L(A0)^/L^ ~ 1/L, thus 
vanishing in the thermodynamic limit. Thus, the isotropic-nematic interface with an oblique 
anchoring constraint imposed on the nematic side can be regarded as being marginally 
stable, as opposed to unstable, provided the thermodynamic limit is taken with care. We 
demonstrate that suitably chosen, fiexible variational choices for the uniaxial and biaxial 
profiles can capture the variation of components of the Q tensor as a function of space. 
These variational profiles are obtained by generalizing results from a calculation of biaxial 
and uniaxial order parameter profiles in the planar case. These profiles are benchmarked 
against numerical calculations. 

The outline of this paper is the following. In Section [Tll we briefiy review aspects of the 
Landau- Ginzburg-de Gennes transition which will be required in our analysis and obtain 
the equations representing the variational minimum of the GLdG free energy, in a basis 
adapted to the symmetry of the problem. Section IIIII describes solutions to these equations, 
as appropriate to the cases of planar and homoeotropic anchoring. The classic tanh profile 
obtained by de Gennes is an exact representation of the interface in the limit of homoeotropic 
anchoring as well as when the L2 elastic constant vanishes, in which case the interface is 
stable for any anchoring condition. In Section II VI we present our numerical approach to the 
problem of interface structure, showing how numerically exact profiles for the variation of 
S, T and 6 can be obtained within the framework of a minimization of the full GLdG free 
energy, subject only to the condition that an interface is forced into the system. 

In Section lYl we describe our variational approach to this problem, motivating the choice 
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of a three-parameter variational ansatz inspired by the approximate solution due to Popa- 
Nita, Sluckin and Wheeler. We show that this variational ansatz captures the features 
of the solution in both the extreme cases of planar and homoeotropic anchoring, and is 
flexible enough to describe the intermediate regime as well. In Section Wl\ we describe 
our methods of minimization for the variational problem and our results for L2 > and 
L2 < 0. We describe how our numerical and variational calculations can be used to provide 
an accurate picture of the interface with an oblique anchoring condition In Section IVIII we 
present asymptotic results for the variation of S, T and 6 close to the bulk nematic state. 
Section IVIIII contains our conclusions. 



II. THE GINZBURG-LANDAU-DE GENNES APPROACH TO THE ISOTROPIC- 
NEMATIC TRANSITION 



The Ginzburg-Landau-de Gennes free energy functional F = Fh + F^i [5| is obtained from 
a local expansion in powers of rotationally invariant combinations of the order parameter 
Q(x,t), 

J-,[Q] = ^ATrQ' + ^BTrQ' + ^C{TrQ'f + E'iTrQ'f (2) 

The restriction to the terms shown above are sufficient to yield a first-order transition 
between isotropic and nematic phases as well as a stable biaxial phase, obtained when 
E' ^ o|4|. 

To this local free energy, non-local terms arising from rotationally invariant combinations 
of gradients of the order parameter must be added. The choice of the following two lowest- 
order gradient terms is common 

^el[dQ,] = ^Li{daQf3'y){daQf^'y) + ^L2{daQ al3){d^Q p-y) , (3) 

where 7 denote the Cartesian directions in the local frame, and Li and L2 represent 
the elastic cost for distortions in Q[4]. The fact that there are only two terms which appear 
to this order implies that only two of the three Frank constants are independent. The limit 
in which L2 = 0, or of zero elastic anisotropy corresponds to the case in which all Frank 
constants are equal. The relationship between Li and L2 and the Frank constants Ki,K2 
and K3 are the following: Ki = K2 = 9/4(2Li + Ls)^^ ^nd K3 = 9/2LiS^^, 4\. Note that 
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K = L2/L1 negative is allowed, although k < 1.5 must be satisfied to ensure positivity of 
the elastic constants. 

In the free energy density of Eq. [21 A = Aq{1 — T/T*), where T* denotes the supercooling 
transition temperature. From the inequality i(TrQ^)'^ > (TrQ^)^, higher powers of TrQ^ 
can be excluded for the description of the uniaxial phase. Thus the uniaxial case is described 
by £" = whereas E' ^ for the biaxial phase. We will assume that E' = 0, thus ensuring 
that the stable ordered phase is the uniaxial nematic. For nematic rod-like molecules B < 
whereas for disc-like molecules, B > 0; for concreteness, we will assume B < here. The 
quantity C must be positive to ensure stability and boundedness of the free energy in both 
the isotropic and nematic phases. 

The first order isotropic to uniaxial nematic transition at the critical value = is thus 
obtained from, 

A = ^CS! (4) 
B = -^CS,. (5) 

We choose B = —0.5, C = 2.67 and A = B'^/27C, thus enforcing phase coexistence between 

n 

an isotropic and uniaxial nematic phase 

The interface is taken to be fiat and infinitely extended in the x—y plane. The spatial vari- 
ation of the order parameter only occurs along the z direction j^. We scale Qa/B ~^ Qap/Sc 
where Sc = ^ g^^*-^' ^^"^ measure lengths in units of Ic = a/54C(Li + 2L2/3)/i?2. 



A. The Ginzburg-Landau-de Gennes Equations 

The director n, sub-director 1 and their joint normal m together define a frame. We 
define z as the direction perpendicular to the interface. The fixed orientation of the nematic 
axis at z ^ 00 can be used to define a plane, the xz plane. From symmetry, and following 
the arguments of Sen and Sullivan, the nematic director must always remain in this plane j^. 
Thus, the spatial dependence of the frame orientation can only come from the variation of 
a single tilt angle 6', which is measured between the z axis and n, as shown in Fig. [H 

Since we assume a fiat interface, the components of Q are functions only of z. The tensor 
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Q n the local frame defined by the principal axes, is diagonal and given by 



Q 



V 



-{S + T)/2 
-{S-T)/2 

5 



(6) 



Transforming to the space-fixed frame (the laboratory frame), by rotation through the ap- 
propriate angle yields 



Q6 



v 



COS 6 sm.6 

1 
- sin ^ cos 



\ 



( 



V 



-{S + T)/2 

-{S-T)/2 

5 



\ 



V 



cos 6 — sin^ 

1 
sin^ cos^ 



(7) 



■i(35 + T) sin 2^ 




(8) 



Thus, Qt) takes the form 

( -\{S + T)cos^e + Ssm^e 

-{S-T)/2 
-\{3S + T)sm2e + T) sin^ ^ + -^cos^ ^ J 

Inserting this tensor form into the elastic free energy ^ei[Q] yields the elastic contribution 
to the free energy 

(12 + 5A; + 3kcos(2e))d,S'^ + Aksm^{d)d,Sd,T + 2 (2 + ks\u^{e)) d,T^ 



8(3 + 2A;) 



2A;sin(2^)(35 + T) {d,S - d^T) d^9 (2 + k){3S + Tfd.O^ 



+ 



(9) 



8(3 + 2A;) ' 4(3 + 2A;) 

Note that this contribution must vanish if T and 6 are uniform. 

The bulk free energy contribution ^/i[Q] is unchanged, as a consequence of the fact that 
the Landau term is constructed from rotationally invariant terms in the order parameter. It 
then takes the form 

3 9 
The Euler-Lagrange equations minimizing the full free energy F, are obtained from 

_ aF> ^ aF> 

where 9 = dO/dz. This yields 



4(2 + k) {ZS' + T') 9' -ksm{29){S" - T") + 2(2 + k){ZS + T)9" = 0, 



(12) 



which can further be simphfied as 

- k sm{20){3S + T){S" - T") + (2(2 + k){?,S + TfO')' = 0, (13) 

where the primes indicate derivatives with respect to z. 

First, note that for k = {i.e. no elastic anisotropy) the above equation has only the 
solution 6' = 0, implying that 6 is constant. A similar situation holds for the special 6 values 
^ = 0, 90°, for which again the only solution has 6' = 0. Thus, in these special limits, the 
angle 6 remains fixed throughout the system. These results are, of course, consistent with 
the result that planar {6 = 90°) and homoeotropic {6 = 0) anchoring conditions yield a 
well-defined interface. Also, provided elastic anisotropy is absent, one can continue to define 
a stable interface for an arbitrary 6, since 6 sticks to its asymptotic value throughout. 

Finally, we note that once S and T are saturated, S' = T' = S" = T" = 0, and thus 6' = 
constant, yielding a linear variation of 6 with z. 

For completeness, the full set of Euler-Lagrange equations representing the minimization 
of the GLdG equations are, in addition to the 6 equation above 

_ nkcos{2e) + 6 + 3k)i3S + T)\ ^ ^ k{4 + 3k + k cos{2e)) sm{ef 
~\ 6 + 4k ) 4 (6 + 7A; + 2k^) 

96 + 88fc + 19k' + 12fc(2 + k) cos(2g) + P cos(4g) \ ^„ ^ _ ^ ^ ^ ^^T^ 
16 (Q ~h Tk ~\- J 3 

(14) 





32 + 24k + 3^ - 4k{2 + k) cos{29) + P cos(4e) ^^ ^„ ^k{4 + 3k + k cos(2^)) sin^OfS" 



16 (6 + 7A; + 2A;2) J 4 (6 + 7fc + 2fc2) 

, {kcos{2e)-2-k){3S + T)e'' _2^ , , 4 . , 4^,^ 



-T + 4ST + -T' + -S^T 



6 + 4k 3 9 3 

(15) 



III. INTERFACE STRUCTURE FOR PLANAR AND HOMOEOTROPIC AN- 
CHORING 

This section briefly reviews the methodology for the solution of interfacial structure in 
the cases of homoeotropic and planar alignment . While the exact solution in the case 
of homoetropic alignment, as originally proposed by de Gennes, motivates the canonical 
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tanh form for the uniaxial order parameter, the more complex situation of planar anchoring 
requires the simultaneous solution of equations of motion for both S and T, in addition to 
the equation for 9 6|. We discuss how the Popa-Nita, Sluckin and Wheeler solutionj^ of the 
planar case can be generalized, in a variational sense, to the more general problem of an 
oblique anchoring condition. 

A. Homeotropic Alignment 

The equation of motion for homoeotropic boundary conditions is easily obtained by set- 
ting 6' = 0, in the defining equations above. This immediately yields, 

-d^S = S-3S^ + 2S'' + T' + (16) 
2 3 

dlT = tT + 2ST+^- + (17) 



2(3 + 2fc) ^3 9 3 

It is easy to see that these equations have the solutions 



5 = l(l + tanh(-|)), T = 0; (18) 
Here the treatment of de Gennes is exact. 

B. Planar Alignment 

The case of planar alignment follows from setting 9 = 90° in the Euler-Lagrange equa- 
tions. This then yields the following set of coupled partial differential equations for the S 
and T order parameters. 



(3 + 2fc) " (3 + 2fc) " 3 

* .9js + ^c).^r = ir + 8Sr + ?^! + ?^. (20) 



(3 + 2fc) ^ (3 + 2A;) ^ 3 9 3 

In the zeroth order aproximation we drop terms in T as in the solution of the first equation. 
This then yields 5* = ^(1 + tanh(^^)) where ^ = -y/^^^- Putting this in equation (l20j) . 
scaling z again with \/2^ and neglecting the nonlinear term, we get the following equation. 

d^T = 2p{tanh^{z) + 8tanh(2) + 9)T 
k 

H -tanh(z)(l + tanh(z))(l -tanh(2)), (21) 
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with /? = 3^1). 

The PSW approximation now consists of dropping the d'^T term, yielding the algebraic 
equation 

k 

2(3{tanh'^{z) +8tanh(z) + 9)T = tanh(2)(l + tanh(2))(l - tanh(z)), (22) 

2 + k 

which then immediately yields 

^ k tanh(2;)(l + tanh(2;))(l — tanh(2;)) . . 

~ ~2(3{2 + k) (tanh2(z) +8tanh(2) +9) ' ^ ^ 

We have recently suggested an improvement to these results, motivated by our tests of 
the self-consistency of the PSW approximations. These tests indicate that the d'^T term 
dropped by PWS should be retained for a more accurate description of the interface. Our 
analytic results for this case, expressed as a sum over hypergeometric functions, agree well 
with numerical solutions of the GLdG equations and represent a significant improvement 
over the PSW solution, particularly in the case of small k. 

IV. NUMERICAL MINIMIZATION OF THE GINZBURG-LANDAU-DE 
GENNES FREE ENERGY FOR THE INTERFACE PROBLEM 

Our numerical results for the isotropic-nematic interface with an oblique anchoring condi- 
tion are obtained from a direct minimization of the Ginzburg-Landau-de Gennes functional, 
with boundary conditions which ensure the presence of the interface as well as impose the 
required anchoring condition on the 6 field. Our numerical methodology is the following: 
Defining a system size L, we discretize the one-dimensional (z) coordinate into points, 
defining 6 = L/N. We use, typically, N = 1001. The values of the fields S, T and 9 at each 
of these points is varied so as to minimize the combined integrals of Eq. [9] and Eq. [TOl 

To do this, we perform a straightforward evaluation of the integral using the trapezoidal 
rule, replacing derivative terms in the integrand by the finite difference approximants. Thus, 
the gradient term dS/dz ~ [S{i + 1) — S{i)] /6. We have also used a variable discretization 
in some of our calculations, to assess the accuracy of our results, sampling with closely spaced 
points in the vicinity of the interface where the variation of S and T is largest. We impose 
boundary conditions on S, T and 6, by fixing the values at the two extreme boundaries to 
their values in the isotropic (S* = 0, T = 0) limit, with 6 arbitrary, and in the nematic limit 

{s = i,T = o,e = e,). 
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FIG. 2: [Color online] Profiles of the biaxial (T) (main figure) and uniaxial (S) order (inset) 
parameter as a function of the coordinate z across the interface, for planar anchoring and k = 4 
as obtained from a direct numerical minimization of the LGdG functional (solid line). The results 
obtained from an spectral collocation method are shown as points. 

The location of the interface is fixed at the centre, by imposing 5* = 1/2 at the central 
site. In principle, in a system of finite size L, our methods yield a constrained minimum 
for the following reason: The elastic energy on the nematic side is minimized by allowing 
the nematic region to expand as far as possible, effectively forcing the interface to invade 
the isotropic side. However, as discussed above, in the thermodynamic limit L —>■ oo, this 
elastic energy cost reduces as 1/L, vanishing in the thermodynamic limit where a stable 
interface is obtained. Alternatively, one can think of this in terms of adding a localized 
pinning potential with strength vanishing as L —>■ oo, which serves only to stabilize the 
location of the interface. 

This relatively simple approach yields results of very high quality, as we have checked by 
a direct comparison to exact results for the planar anchoring well as to numerical 
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calculations using spectral methods in the case of planar anchoring. We have used the 
minimization routines (NMinimize) in Mathematica to find the stationary values of S,T 
and 6 which minimize the free energy subject to the applied boundary conditions. This 
routine selects the most appropriate methodology from a variety of minimization techniques 
available, iterating till an accuracy between successive iterations of 1 part in 10*^ is obtained. 

As a test of the quality of the minimization methodology which will be used in this 
paper, we show in Fig. [2l profiles of the biaxial (T) (main figure) and uniaxial {S) order 
(inset) parameter as a function of the coordinate z across the interface, as computed by 



the numerical spectral methodology of Ref. 22| and the minimization technique described 



above, for the case of planar anchoring i.e. 9^ = 90°, with k = 4. Results obtained from the 
numerical minimization of the LGdG functional are shown as the solid line whereas results 
from the spectral collocation scheme of Ref. 



22| are shown as points. These coincide to high 



accuracy. 



V. VARIATIONAL METHOD 



Clearly, the solution of the full set of equations for 5*, T and 6 given above is a formidable 
problem. Our approach to this problem therefore proceeds through the construction of 
simple, physically motivated variational choices for 6{z), S{z) and T{z). This choice is made 
keeping in mind that requirement that the results should be consistent with computations 
in the simpler ^ = 0, 90° limits, where the angular variation is absent and the de Gennes 
solution and the PSW solution are obtained, respectively. 

Our approach begins by assuming a profile of the form 



w NN 1 , ,/ , ( 1 + t auh (cz) ) ( 1 — t auh (c2;) ) 

-(1 + tanh(az)) and T = -btanh(cz)- ^ 

2' tanh^(c2) + 8tanh(c2) + 9 



together with the assumption that the theta variation can be fitted to a simply parametriz- 
able function. We have examined a variety of such functions for the case of planar anchoring, 
including (a) ^ = 90° - 2^z for ^ > 0, 90° for z < 0, (h) 6 = 90° - f (1 + tanh(aiz)) which 
implies that at ;z = oo, 6* = 90° - and at ;z = -oo, 9 = 90°, (c) 6* = 90° - f (1 + tanh(ai;z)) 
which implies that at z = oo, 6 = 90° — ip and at z = — oo, 6 = 90°, (d) 6 = ^(1 + tanh(aiz)) 
(e) e = ^z + iIj/2 and (f) ^ = p + |(1 + tanh(ai;z)). 
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Our best results are obtained with the variational form 

= p + ^py (25) 

subject to a constraint p + ip = 6l where 6l is the value of angle at L, the system size. It 
will be our intention to take the L oo limit later. 

Note that the choice p = 90°, = 0, a = l,b = 2p{2+k) I'scovers the profile of PSW for 
the planar case. The parameter values ip = 0,b = generate the de Gennes solution. Thus, 
the two extreme limits of the variation of the anchoring angle can be obtained with the 
appropriate choice of parameter values in the variational form chosen above. These can be 
simply generalized to the case of homoeotropic anchoring. 

VI. NUMERICAL METHODOLOGY FOR THE VARIATIONAL SOLUTION 

These variational ansdtze for S and T are inserted into the form for the free energy, 
which is then minimized with respect to the parameters a, b, c, and p. This minimization 
is carried out using Mathematica. We use the " Nelder-Mead" method for the minimization 
of a function of n variables. This is a direct search method which uses an initial choice of 
n + 1 vectors which form the vertices of a polytope in n— dimensions and a methodology for 
changing the vertices of this polytope iteratively. The process is assumed to have converged 
if the difference between the best function values in the new and old polytope, as well as 
the distance between the new best point and the old best point, are less than preset values, 
typically of the order of 10~^°. 

To eliminate problems arising from an incorrect choice of initial values, we have computed 
the minima for about 100 separate initial conditions and chosen the parameter values corre- 
sponding to the least value of the free energy from these. Our results for the minimization 
have been crosschecked using the differential evolution method, a simple stochastic function 
minimizer. 

A. Results from the Numerical and Variational Minimization: k > 

In Fig. [3], we show profiles of the biaxial (T) and uniaxial (S) order parameter as a 
function of the coordinate z across the interface. These are computed by direct numerical 
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minimization of the LGdG functional, via the methodology described in the previous section. 
We allowed 6 on the isotropic side to vary, finding that the free energy minimum was obtained 
when 6 was stuck to the value it attained at the location of the interface. This value is 
somewhat smaller than 90° for small system sizes but asymptotes to this value as L goes to 
infinity. 

We show the T profile in the main sub-figure for systems of size L = 50, 1290 and 
parameter values (a) k = 1,6^ = 30 (b) k = 8,6e = 30, (c) k = l,6e = QO and (d) 
K = 8, 6'e = 60. N and I in the figure refer to nematic and isotropic respectively. The insets 
to each of (a), (b), (c) and (d) show the corresponding profiles for S. 

We note that for larger anchoring angles, the T profile converges faster as a function of 
system size than for smaller angles; contrast the behavior for 6^ = 30° and 6^ = QO" in the 
figure. The profiles are qualitatively similar to profiles obtained for the 6^ = 90° degree, and 
asymptotically match this profile as L oo. 

In the inset to Fig. HJ we show the profile of 6, the angle describing the orientation 
of the local director field as a function of z across the interface, as obtained from our 
numerical minimization. We show data for systems of size L = 125,258,516 and 1290, and 
for parameter values (a) k = 1,9^ = 30 (b) k = 8,9f, = 30, (c) k = 1, 6'e = 60 and (d) 
K, = 8,6e = 60. The main figure, in each case, plots the same data as a function of the 
scaled coordinate z/L on the x— axis and the quantity {9 — 9e)/{9,nax — ^e) on the y— axis , 
thus normalizing the value to its maximum. This produces high quality collapse of the data, 
indicating that the angle profile is linear on the nematic side, interpolating linearly between 
its value at the interface to the anchored value of 9^ a.t z = L. Also, as the system size is 
increased, the value at the interface {z = 0), approaches 90°, indicating that anchoring at 
the interface is always planar in the asymptotic limit. 

In Fig. [5] we show the comparison between the computed 3-parameter variational profile 
for the angle ^ as a function of the coordinate z across the interface, for a system of size 
L = 125, as obtained from a direct numerical minimization of the LGdG functional (solid 
line) and from the variational calculation described in the text (point). These are shown 
for parameter values (a) k = 1,9^ = 30 (b) k = 8,9^ = 30, (c) k = 1,9^ = 60 and (d) 
K = 8, = 60. The inset labeled (i) in each sub- figure shows the corresponding profile of 
S, whereas the inset labeled (ii) shows the profile of T. Note that the variational result 
coincides with the result obtained from a direct numerical minimization to high accuracy. 
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FIG. 3: [Color online] Profiles of the biaxial (T) and uniaxial (5) order parameter as a function 
of the coordinate z across the interface, computed by direct numerical minimization of the LGdG 
functional. These are shown in the main figure for systems of size L = 50, 1290 and parameter 
values (a) k = 1, 6*6 = 30 (b) K = 8,ee = 30, (c) k = 1, = 60 and (d) K = 8,de = 60. The insets 
to each of (a), (b), (c) and (d) show the corresponding profiles for S. N and I refer to nematic and 
isotropic respectively. 
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FIG. 4: [Color online] Main figure: Scaled profiles of the angle 9 describing the orientation of the 
local director field as a function of the scaled coordinate z/L across the interface, as obtained from 
direct numerical minimization of the LGdG functional for systems of size L = 125, 258, 516 and 
1290. These are shown for parameter values (a) k = 1, = 30 (b) k = 8,0e = 30, (c) k = l,6e = 60 
and (d) k = 8,8e = 60. The insets to each of (a), (b), (c) and (d) show the corresponding unsealed 
profiles for 6. 
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As the system size is increased, the value of 9 at the interface approaches 90° within both 
the variational and the direct numerical minimization approaches. 

Fig. O shows the variational parameters a (a), h (b) and c (c) as a function of system 
size L, together with the variation of the variational angle p (d), plotted for k = 1. These 
parameters converge to their L ^ oc values corresponding to the case of planar anchoring. 
In all cases the parameter p converges to the asymptotic value of 90° as the system size is 
increased, consistent with planar anchoring. 

B. Results from the Numerical and Variational Minimization: k < 

Stability imposes the requirement that 3 + 2k > 0, but does not constrain the sign of k (or, 
equivalently L2), apart from this requirement. In this section we explore the consequences 
of a negative value for L2. 

We find that, consistent with de Gennes' prediction, a negative L2 ( or k) consistent with 
stability favours homoeotropic anchoring at the interface, in contrast to the case of positive 
L2. Thus, the biaxiality T generically vanishes as L ^ 00, whereas S assumes the canonical 
tanh form obtained by de Gennes. This can be seen from Fig. [7] which shows the variation 
of S and T, for L = 125, plotted for k, = —1. The anchoring at L is set to an oblique angle 
of 30°. The S and T profiles are consistent with T = for homoeotropic anchoring. 

The preference for homoeotropic anchoring can be seen from Fig. [8] which shows the 
director tilt angle scaled to its minimum value for each system size {L = 125, 258 and 516, 
against z/L for k = —1, where an asymptotic, oblique anchoring angle of 30° is imposed 
on the system at L. The inset shows the bare angles as a function of z for these different 
system sizes. The excellent data collapse indicates that angle profiles in the case of L2 < 
scale in the same way as the L2 > case, except that homoeotropic anchoring is favoured 
in this case. 

Finally, in Fig.[9|, we show, in the main figure, the profile of T, the biaxial order parameter, 
for K = — 1, in the two extreme cases of planar (0°) and homoeotropic (90°) anchoring, with 
L = 50. Importantly, the profile of T is inverted with respect to profiles obtained for k > 0, 
with the minimum appearing on the isotropic side of the interface rather than the nematic 
side, as earlier. The profile of S is consistent with a tanh form. While the profile of T is 
non-zero for planar anchoring, biaxiality vanishes for the homoeotropic anchoring case. 
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FIG. 5: [Color online] Main figure: Profiles of the angle 6 describing the orientation of the local 
director field as a function of the coordinate z across the interface for a system of size L = 125, 
as obtained from a direct numerical minimization of the LGdG functional (solid line) and from 
the variational calculation described in the text (point). These are shown for parameter values (a) 
K = 1,6*6 = 30 (b) K = 8,0e = 30, (c) K = 1,6*6 = 60 and (d) k = 8,6*6 = 60. The inset labeled 
(i) in each sub-figure shows the corresponding profile of S, whereas the inset labeled (ii) shows the 
profile of T. N and I refer to nematic and isotropic, respectively. 
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FIG. 6: [Color online] The variation of the variational parameters a (a), 6 (b) and c (c) with system 
size L, together with the variation of the variational angle p (d), plotted for k = 1. Note that these 
parameters quickly converge to their L ^ oo values corresponding to the case of planar anchoring. 
In all cases the parameter p appears to converge to the asymptotic value of 90° as the system size 
is increased. 
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FIG. 7: [Color online] The variation of the S and T, for system size L = 125, plotted for k = — 1, 
with an oblique anchoring angle of 30°. Our results are consistent with T = for homoeotropic 
anchoring. 

These results are consistent with the general trends observed in the case of k > 0, with 
the difference that homoeotropic, rather than planar, anchoring is preferred once k turns 
negative. 



VII. ASYMPTOTIC SOLUTION 



We can use our ansatz for S and T to check the self-consistency of our conjectured 
behaviour for 6. Our chosen forms imply S = 1 — e'"^""^ and T = — 6e~^"^ deep into the 
nematic phase, as z ^ oo. Then 5" = 2ae~^"^, T' = 2abe~'^°'^ and 5"' = — 4a^e~^"^, 
T" = — 4a^6e~^"^. Inserting these into the equation for 9 as below, 

4(2 + k) (35' + T) 6' - A;sin(2^)(5" - T"), +2(2 + k){?,S + T)e" = 0, (26) 
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FIG. 8: [Color online] Plot of the angle scaled to its minimum value for each system size (L = 
125,258 and 516, against z/L for k = —1 and an asymptotic, oblique anchoring angle of 30°. 
The inset shows the bare angles as a function of z for these different system sizes. Note that the 
excellent data collapse indicates that angle profiles in the case of L2 < scale in the same way as 
the L2 > case, with a homoeotropic anchoring being favoured at the interface. 

we get 

8(2 + A;) (3 + h) ae^^^'e' + k sin(2^) (1 - h)a^e-^''' + 6(2 + k)e" = 0. (27) 

As z — s> 00, this equation reduces to 9" = 0. Thus, 9 should have a linear profile in this 
asymptotic limit, taking the form 

9 = p + iJY. (28) 

We can also compute corrections to this profile for z — >■ oo~. Let us now expand about 
the z = 00 limit, in which case 9" = 0. Thus, 

9' -ka{l-b) 



sm{29) 2(2 + A;) (3 + 6)' 
Integrating the left-hand side of this equation, we obtain 



(29) 



l|ntanW-lnC^ ^-/;';^-^>;^ . (30) 
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FIG. 9: [Color online] Main Figure: Profile of T, the biaxial order parameter, for k = —1, in the 
two extreme cases of planar (0°) and homoeotropic (90°) anchoring. Note that the profile of T is 
inverted with respect to profiles obtained for k > 0, with the minimum appearing on the isotropic 
side of the interface. INset: The profile of S, consistent with a tanh form. Data are computed 
for L = 50. While the profile of T is non-zero for planar anchoring, biaxiality vanishes for the 
homoeotropic anchoring case. 

— ka(l — b)z 

which has a solution 6 = tan~^[Ce (2+^)13+6)]. It can be seen that this will vanish as z goes to 
oo and is, in effect, negligible apart from a region close to the interface, at z = 0. 

VIII. SUMMARY AND CONCLUSIONS 

In this paper, we have presented our results for the problem of the isotropic-nematic inter- 
face within Ginzburg-Landau-de Gennes theory, for the case in which an oblique anchoring 
condition is imposed on the system asymptotically on the nematic side. In this case, we find 
that nematic elasticity dictates that the nematic orientation smoothly interpolates between 
a value of 90° at the interface (planar anchoring) to the anchored value at the boundary on 
the nematic side when k > 0. Thus, the preferred value of the anchoring angle at the inter- 




24 



face is 90° in this case. The case k, <0 with k, satisfying the stabihty requirement k, > —1.5 
leads to stable homoeotropic anchoring at the interface, as predicted by de Gennes. 

We have used simple variationally based descriptions of the structure of the interface, 
with our methods capturing essential features of interface structure, both qualitatively and 
quantitatively, for the case of oblique anchoring. Our methods access the non-trivial struc- 
ture of biaxiality at the interface, including the large tail towards the isotropic side and the 
change in the sign of the biaxial order parameter across the interface. Our approach also 
captures the inversion of the profile of biaxiality as k crosses zero. 

The results presented here are broadly consistent with results from density functional 
approaches, molecular simulations and approaches based on the Onsager functional, but 
necessitate fewer approximations, truncations or assumptions about specific model systems. 
Thus, coarse-grained approaches based on the Ginzburg-Landau-de Gennes functional pro- 
vide a powerful methodology for understanding generic features of the isotropic-nematic 
interface. 
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